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Abstract 


This  report  summarizes  our  progress  and  accomplishments  for  the  project  “Secure 
Cooperative  Communications  and  Testbed  Development”,  funded  by  AFRL  under  project 
FA8750-06-2-0167  from  May  8,  2006  to  January  8,  2007.  We  have  worked  on  physical-layer 
security  and  cooperative  communications,  and  have  been  developing  a  wireless  cooperative 
transmission  testbed  to  demonstrate  our  theoretic  results.  We  have  finished  or  made  significant 
progress  on  all  topics  from  the  original  proposal. 

In  this  project,  we  have  obtained  some  major  achievements.  We  have  developed  an  innovative 
secure  transmission  scheme,  have  studied  the  limit  of  cooperative  transmissions,  and  have 
invented  a  new  cooperative  OFDM  transmission  scheme  to  combat  transmission  asynchronism. 
They  are  helpful  to  the  development  of  future  physical-layer  wireless  information  assurance 
techniques  as  well  as  the  cooperative  communication  techniques.  We  have  successfully 
implemented  the  wireless  cooperative  transmission  testbed,  implemented  array  transmissions, 
and  estimated  channels  which  verified  partially  the  validness  of  our  theories. 

First,  as  shown  in  Part  I,  for  secure  transmissions,  we  have  proposed  the  use  of  signal 
processing  techniques  to  protect  wireless  transmissions  as  a  way  to  secure  wireless  networks  at 
the  physical  layer.  This  approach  addresses  a  unique  weakness  of  wireless  networks  whereby 
network  traffic  traverses  a  public  wireless  medium  making  traditional  boundary  controls 
ineffective.  Specifically,  a  randomized  array  transmission  scheme  is  developed  to  guarantee 
wireless  transmissions  with  inherent  low-probability-of-interception  (LPI).  In  contrast  to 
conventional  spread  spectrum  or  data  encryption  techniques,  this  new  method  exploits  the 
redundancy  of  transmit  antenna  arrays  for  deliberate  signal  randomization  which,  when 
combined  with  channel  diversity,  effectively  randomizes  the  eavesdropper’s  signals  but  not  the 
authorized  receiver's  signals.  The  LPI  of  this  transmission  scheme  is  analyzed  via  proving  the 
indeterminacy  of  the  eavesdropper's  blind  deconvolution.  The  proposed  method  is  useful  for 
securing  wireless  transmissions,  or  for  supporting  upper-layer  key  management  protocols. 

Second,  as  shown  in  Part  II,  for  cooperative  transmissions,  in  contrast  to  most  existing  work 
that  depends  on  the  assumption  of  perfect  synchronization  among  the  transmitters,  we  consider 
the  inevitable  asynchronism  among  the  cooperative  transmitters.  Our  studies  show  that  the 
interference  caused  by  asynchronism-induced  channel  dispersion  introduces  an  average  signal- 
to-interference-ratio  (SIR)  degradation  of  at  least  10  dB  in  flat  fading  environment,  which  may 
greatly  limit  the  cooperative  diversity  gain.  We  also  show  that  in  OFDM  transmissions,  the 
cyclic  prefix  (CP)  can  be  used  to  resolve  such  asynchronism.  In  particular,  for  the  first  time,  we 
show  that  the  carrier-frequency  offset  (CFO)  can  be  removed  completely  with  extended  CP.  This 
finding  is  very  useful  for  cooperative  OFDM  transmissions  and  for  the  future  multi-carrier-based 
OFDM  transmissions. 

Finally,  in  Part  III,  a  testbed  is  setup  in  order  to  demonstrate  the  secure  transmission  scheme 
and  the  cooperative  transmission  scheme,  from  which  many  important  assumptions  have  been 
demonstrated,  such  as  transmission  asynchronism,  channel  dispersion,  channel  independence,  etc. 
We  implement  BPSK/QPSK  array  transmissions  using  ComBlock  modules.  Channels  are 
measured  and  used  in  the  simulation  of  secure  transmissions. 
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This  research  has  generated  4  journal  papers  (one  formally  accepted,  one  in  2nd  round  review, 
two  in  preparation),  and  4  conference  papers  published  in  major  conferences. 

In  summary,  we  have  made  significant  progress  and  achievements  in  this  project,  and  thus 
satisfactorily  conclude  this  project.  Nevertheless,  there  are  still  many  interesting  problems  open 
which  deserve  continuous  investigation  should  future  funding  be  available. 
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Background 


This  is  the  final  report  for  the  project  “Secure  cooperative  communications  and  testbed 
development,”  funded  by  Air  Force  Research  Laboratory  under  grant  FA8750-06-2-0167. 

This  report  consists  of  three  parts.  Part  I  develops  physical-layer  security  transmission 
techniques.  Part  II  addresses  cooperative  communications.  Part  III  is  devoted  to  the  testbed 
development.  To  save  space,  many  details  have  to  be  skipped,  but  can  be  referred  from  our 
publications  listed  below: 

[PI]  X.  Li,  J.  Hwu  and  E.  P.  Ratazzi,  “Using  antenna  array  redundancy  and  channel 
diversity  for  secure  wireless  transmissions,”  submitted  to  Journal  of  Communications 
(formally  accepted). 

[P2]  X.  Li  and  F.  Ng,  “Carrier  frequency  offset  mitigation  in  asynchronous  cooperative 
OFDM  transmissions,”  in  2nd  round  review  in  IEEE  Trans.  Signal  Processing. 

[P3]  F.  Ng  and  X.  Li,  “CFO-resistant  receiver  for  asynchronous  MC-DS-CDMA 
systems,”  Proceedings  ofICASSP'2007,  Honolulu,  Hawaii,  April  15,  2007. 

[P4]  X.  Li  and  F.  Ng,  “Using  cyclic  prefix  to  mitigate  carrier  frequency  and  timing 
asynchronism  in  cooperative  OFDM  transmissions,”  in  Proceedings  of  the  40th  Asilomar 
Conference  on  Signals,  Systems  and  Computers ,  Pacific  Grove,  CA,  Oct.  29-Nov.  1, 

2006. 

[P5]  X.  Li  and  J.  Hwu,  “Performance  of  cooperative  transmissions  in  flat  fading 
environment  with  asynchronous  transmitters,”  IEEE  Military  Communications 
Conference  (MILCOM  2006),  Washington,  DC,  Oct.  23-25,  2006. 

[P6]  X.  Li,  J.  Hwu  and  E.  P.  Ratazzi,  “Array  redundancy  and  diversity  for  wireless 
transmissions  with  low  probability  of  interception,”  Proceedings  of  ICASSP'2006, 
Toulouse,  France,  May  14,  2006. 
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Part  I:  Using  Antenna  Array  Redundancy  and  Channel 
Diversity  for  Secure  Wireless  Transmissions 


1.  Introduction 

Along  with  the  rapid  development  of  wideband  wireless  communication  networks,  wireless 
security  has  become  a  critical  concern.  Compared  with  wireline  networks,  wireless  networks  lack 
a  physical  boundary  due  to  the  broadcasting  nature  of  wireless  transmissions.  Any  receivers 
nearby  can  hear  the  transmissions,  and  can  potentially  record  or  analyze  the  transmitted  signals, 
or  conduct  jamming.  This  makes  wireless  security  design  a  challenging  task,  and  the  challenge 
becomes  even  more  severe  if  considering  together  other  unique  characteristics  of  wireless 
networks,  such  as  severe  energy/bandwidth  constraints  of  wireless  nodes,  unreliable/untrusted 
wireless  links,  and  dynamic  wireless  network  topology.  Noticing  that  the  challenge  is  closely 
related  to  the  unique  physical-layer  of  wireless  communications,  physical-layer  security 
techniques  are  thus  helpful,  since  they  can  be  more  effective  in  resolving  the  boundary, 
efficiency,  and  link  reliability  issues. 

One  of  the  important  objectives  of  physical-layer  security  design  is  to  guarantee  wireless 
transmissions  with  low-probability-of-intercept  (LPI).  In  particular,  we  are  interested  in  LPI 
techniques  which  do  not  directly  rely  on  upper-layer  data  encryption  or  secret  keys. 

Existing  physical-layer  LPI  techniques  can  be  classified  into  three  categories:  i)  Signal  power 
approaches  like  beamforming  and  directional  transmissions,  ii)  scrambling  code  approaches  like 
spread-spectrum,  and  iii)  propagation  channel  approaches.  Traditionally,  spread  spectrum  is  the 
most  widely  used  technique  for  LPI.  However,  when  data  transmissions  are  evolving  toward 
wideband,  spread  spectrum  alone  may  not  be  enough  because  of  the  reduced  space  of  spreading 
gain. 

In  general,  the  security  of  most  existing  approaches  depends  on  some  strong  (and  ideal) 
assumptions,  such  as  eavesdroppers  have  null-receiving  energy,  or  have  no  infonnation  about  the 
spreading  codes,  or  can  not  estimate  the  propagation  channels.  Unfortunately,  these  strong 
assumptions  can  hardly  hold  in  practice.  Beamforming  techniques  can  only  reduce,  but  not 
completely  nullify,  the  signal  energy  toward  eavesdroppers.  Spreading  codes  may  be  easily 
estimated  by  eavesdroppers  from  their  received  signals.  Eavesdroppers  may  use  non-blind  or 
blind  deconvolution  algorithms  to  estimate  channels  and  signals,  which  makes  many  channel- 
based  approaches  such  as  [1]  to  lose  security.  As  a  result,  most  existing  approaches  can  hardly 
guarantee  LPI,  or  can  even  hardly  withstand  a  rigorous  LPI  analysis. 

There  have  been  many  important  advanced  wireless  transmission  techniques  developed  in 
recent  years,  such  as  antenna  array,  channel  diversity  and  channel  deconvolution,  some  of  which 
may  bring  new  opportunities  for  achieving  LPI.  In  [2]-[5],  we  have  shown  that  physical-layer 
security  can  be  realized  based  on  channel  diversity  by  using  antenna  array  transmissions.  This 
idea  in  fact  represents  an  innovative  way  of  secure  waveform  design,  differently  from  the 
conventional  spread  spectrum  or  data  encryption  techniques.  Another  innovative  concept  is  that 
we  rely  on  signal  processing  theory  such  as  the  indeterminacy  of  blind  deconvolution  [6]  for 
security,  rather  than  information  theory  [1][7].  The  advantage  is  that  LPI  can  be  guaranteed 
much  easier  in  more  practical  transmissions. 
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Based  on  [2],  in  this  project  we  propose  a  special  deliberate  randomization  method  for 
designing  the  transmit  antenna  weights.  A  transmission  power  analysis  is  also  conducted  to  guide 
weights  design.  Extensive  simulations  and  the  development  of  a  testbed  are  shown. 

2.  Secure  array  transmission  model 

We  consider  a  wireless  network  where  Alice  transmits  to  Bob  in  face  of  a  passive 
eavesdropper  Eve,  as  shown  in  figure  1.  Alice  uses  J  transmit  antennas  in  the  secure  channel, 
and  may  use  some  other  antennas  communicating  with  Bob  which  form  an  insecure  public 
channel.  This  public  channel  may  be  used  for  the  synchronization  purpose  between  Alice  and 
Bob,  e.g.,  for  Bob  to  track  carrier  frequency  and  timing.  Note  that  such  a  setting  with  a  secure 
channel  and  a  public  channel  is  standard  in  many  information-theoretic  security  studies  or  key 
management  protocols  [7]. 


Figure  1,  Secure  wireless  transmission  model.  Alice  transmits  to  Bob  using  antenna  array, 

in  face  of  passive  eavesdropper  Eve. 


We  consider  only  the  secure  channel  from  Alice  to  Bob.  A  beamforming-like  array 
transmission  procedure  is  used  by  Alice  to  transmit  to  Bob  a  symbol  sequence  {£(/?)}  which  is 
assumed  as  i.i.d.  uniformly  distributed  with  zero-mean  and  unit  variance.  Though  more  complex 
pre-processing  can  be  exploited,  Alice  just  uses  a  simple  weighting  scheme  with  weighting 
coefficients  w,  («) .  Therefore,  through  the  J  antennas,  Alice  transmits  signal  vectors 


s(«)  = 


1 

i _ 

1 

i _ 

Jj(n)_ 

_Wj(n)_ 

b(n)  =  w(/?  )&(/?) 


(1.1) 


where  w^n)  denotes  the  weighting  coefficient  of  the  ;,h  transmit  antenna  during  the  symbol 
interval  n . 


Assume  Rayleigh  flat  fading  channels  and  assume  Bob  use  only  one  receiving  antenna  for 
both  simplicity  and  worst  case  consideration.  Extension  to  receiving  antenna  arrays  can  be  found 
in  [5].  The  signal  received  by  Bob  is 

x(n)  =  'yh’sl(n)  +  v(n)  =  hffs(/?)  +  v(«)  (1-2) 

1=1 

where  v(«)  denotes  AWGN  with  zero-mean  and  variance  crj ,  /?,*  denotes  channel  coefficients 
which  are  independent  complex  circular  symmetric  Gaussian  distributed  with  zero-mean  and  unit 
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variance,  and  h  =  [/?,  /?, ]7 .  (•)*,  (•)'  and  (•)"  denote  conjugation,  transposition  and  Hermitian, 

respectively. 

The  eavesdropper  Eve  may  use  multiple  receiving  antennas  for  better  interception,  and  the 
interception  becomes  much  easier  with  flat-fading  channels.  Therefore,  we  consider  the  worst 
case  to  Alice  and  Bob  where  Eve  receives  signals  from  M  receiving  antennas,  whose  received 
signals  can  be  denoted  as 


x,(«)  =  H,s(«)  +  ve(/?)  (1.3) 

The  vector  ve(/?)  is  AWGN  with  zero-mean  and  covariance  matrix  crlM  ,  where  \M  is  the 
M  xM  identity  matrix. 

We  assume  that  each  element  of  H  has  the  same  distribution  as,  but  is  independent  from, 
those  ofh .  From  the  extensive  studies  on  antenna  array  channels,  we  know  that  as  long  as  the 
distance  between  Bob  and  Eve  is  larger  than  half  of  a  carrier  wavelength,  then  their  channels  can 
be  considered  as  independent  [8]. 

Under  the  above  assumption,  channels  h  and  H  are  different  almost  surely,  especially  when 
J  is  large.  We  further  assume  that  Eve  does  not  know  h  and  H(, .  However,  Eve  may  try  blind  or 
non-blind  methods  to  estimate  H  from  her  received  signal  xe («) .  On  the  other  hand,  Alice  and 
Bob  do  not  know  h  and  H  either.  We  will  discuss  ways  for  Alice  to  obtain  channel  knowledge 
h  since  transmit  beamforming  requires  transmitter-side  channel  infonnation.  Nevertheless,  our 
major  focus  in  this  paper  is  the  design  of  transmission  weights  w(/?)  so  that  Bob  can  detect 
symbols  b(n)  successfully  with  low  bit-error-rate  (BER)  while  Eve  can  estimate  neither  H 
norZ>(«) . 

3.  Array  transmission  with  deliberate  signal  randomization 

To  introduce  high  BER  to  Eve  is  to  prevent  Eve  from  channel/symbol  estimation.  This  means, 
firstly,  Alice  can  not  transmit  training  signals  by  the  J  transmit  antennas,  because  otherwise  Eve 
can  trivially  utilize  such  training  for  channel  estimation  [9]  [10].  Without  training,  the  only  way 
left  for  Eve  is  blind  deconvolution  [6][11][12].  Therefore,  secondly,  Eve’s  blind  deconvolution 
capability  must  be  prevented.  Because  Bob  has  no  more  advantage  over  Eve  on  channel 
estimation,  such  requirements  also  mean  that  Bob  can  hardly  estimate  his  own  channel  h . 

To  meet  both  requirements,  we  propose  a  transmission  scheme  in  which  Bob  can  detect 
symbols  b(n)  without  the  knowledge  of  channel  h  .  In  addition,  we  use  a  deliberate  signal 
randomization  technique  in  this  scheme  to  randomize  Eve’s  signal  but  not  Bob’s  signal  so  that 
blind  deconvolution  of  Eve  has  irresolvable  ambiguity. 

In  order  for  Bob  to  estimate  symbols  b( n ) ,  the  channel  h  from  Alice  to  Bob  has  to  be 
resolved.  In  our  scheme,  we  ask  Alice  instead  of  Bob  to  estimate  and  utilize  the  knowledge  of  h . 
Alice  can  estimate  h  based  on  channel  reciprocity  [8],  where  Bob  first  transmits  a  training  signal 
to  Alice  using  the  same  carrier  frequency  as  the  secure  channel,  from  which  Alice  can  estimate 
the  backward  channel.  Since  the  forward  channel  h  equals  the  backward  channel  according  to 
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reciprocity,  Alice  can  immediately  use  the  estimated  channel  as  h  to  design  transmission 
weights.  Note  that  this  procedure  gives  no  useful  information  to  Eve  because  the  latter  can  only 
estimate  the  channel  from  Bob. 

Our  basic  idea  is  to  make  h"w(«)  a  deterministic  constant,  while  Hw(«)  changing  randomly 
in  each  symbol  interval,  by  exploiting  the  knowledge  of  h .  For  this  purpose,  Alice  designs  the 
transmitting  weights  vector  w(«)  so  that 

hAv(«)=||h||  (1.4) 

Obviously,  if  the  channel  h  is  constant  or  slowly  time-varying,  we  need  J  >2  transmitters, 
which  explains  why  array  transmission  is  necessary. 

From  the  received  signal  x(n)  =||  h  ||  b(n)  +  v(«) ,  Bob  can  detect  symbols  as 


b(n)  =  argmin  |;x(«)—  ||.h  ||  b(n)  \2  (1.5) 

b(n) 

where  ||  h  ||  can  be  easily  estimated  from  the  received  signal  power. 

Alice's  design  of  w (n)  under  the  constraint  (1.4)  can  be  performed  as  follows.  In  each  symbol 
interval  n ,  Alice  first  selects  from  h  randomly  an  element  h  with  sufficiently  large  magnitude. 
The  weighting  vector  w (/?)  is  then  generated  as 


where 


w  (/?)  =  P(«) 


«,-f#"z,(»)' 


«,■  =-r!|h|| 
hj 

«•  - 1  \K  -  V.  -  Kf 

Z  i(n)  =  [w1(n)  •••  ww(n)  wM(n)  •••  w/(n)f 


(1.6) 


(1.7) 


The  matrix  P  («)  is  a  ./  x  j  permutation  matrix  corresponding  to  the  selection  of  ht  from  the 
vector  h,  i.e.,  its  function  is  to  insert  the  first  row  of  the  following  vector  into  the  A  row.  The 
vector  z,0)  is  arbitrary,  whose  dimension  J- 1  is  the  degrees  of  freedom  in  antenna  array 
transmissions  that  we  can  exploit  for  deliberate  signal  randomization. 

This  array  weights  design  procedure  is  outlined  below  as  algorithm  1 . 


Algorithm  1 .  Update  array  weights  vector  w (/?)  in  each  symbol  interval  n 

1 .  Select  randomly  a  channel  coefficient  hi ,  with  sufficiently  large  magnitude  |  /?,  |>  a  . 

2.  Generate  random  variables  Wj(n),  1  <  j<J,j*i. 

3.  Calculate  w(«)  by  (1.6)-(  1.7). 
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The  selection  of  the  threshold  a  in  step  1  will  be  discussed  in  section  4,  whereas  step  2  will 
be  detailed  in  the  sequel. 

One  of  the  major  advantages  of  algorithm  1  is  its  linear  computational  complexity.  Efficient 
computation  is  important  because  w(«)  is  recalculated  in  each  symbol  interval. 

From  (1.6),  we  can  choose  z ,.(/?)  appropriately  to  prevent  Eve  from  blind  deconvolution.  In 
general,  this  purpose  can  be  fulfilled  by  simply  making  z ,.(«)  to  have  a  distribution  unknown  to 
Eve  since  it  is  well  known  that  successful  blind  deconvolution  requires  the  receiver  know  some 
special  statistics  or  structure  of  the  transmitted  signals  [6]  [9].  However,  existing  results  of  blind 
deconvolution  are  mostly  on  how  to  conduct  blind  deconvolution,  not  on  how  to  prevent  blind 
deconvolution.  The  proof  of  the  incapability  of  blind  deconvolution  is  rarely  seen. 

To  furbish  a  rigorous  quantitative  proof  of  the  incapability  of  blind  deconvolution,  we 
consider  a  more  structured  scheme  where  Alice  designs  z ,.(/?)  such  that  r,(«)  =  z,( «)/?(/?)  is  (J-l)- 
variate  Gaussian  distributed  with  mean  n  and  covariance  matrix  £  ,  i.e.,  r,.(«)  □  5Vj-_,(n,  2:) .  The 
parameters  n  and  £  are  arbitrary  and  unknown  to  both  Eve  and  Bob,  and  can  even  be  time- 
varying.  Based  on  this,  a  formal  proof  of  the  LPI  can  be  furbished  by  applying  the  indeterminacy 
of  Eve’s  blind  deconvolution.  Detailed  can  be  found  in  our  paper  [PI].  One  of  the  major 
arguments  is  the  following  proposition. 

Proposition  1.  From  the  distribution  ofx.O) ,  the  channel  matrix  H  is  indistinguishable  from 
HQ  with  a  J  x  j  matrix 


Q  =  P» 


u  0 
0  V 


p ,-'(») 


(1.8) 


where  u  is  an  arbitrary  non-zero  scalar  and  V  is  a  (J-l)x(j-l)  arbitrary  nonsingular  matrix. 


Since  Eve  can  not  estimate  Hc,  what’s  left  for  her  is  a  brute-force  exhaustive  search  of  vector 
h"H;‘  (assume  H  is  invertible).  The  complexity  increases  exponentially  with  the  number  of 
transmit  antennas  J  .  If  Eve  must  use  K  -level  quantization  of  channel  coefficients,  then  the 
brute-force  search  needs  to  consider  at  least  K2J  possible  coefficients  (real  and  imaginary  parts), 
which  means  a  complexity  of 0(K2J).  For  example,  for  QPSK  transmission  at  SNR  25  dB,  in 
order  to  guarantee  BER  0.01,  K  should  be  at  least  128.  In  this  case,  a  J  =  8  transmit  antenna 
array  brings  a  complexity  ofO(2112)  .  This  complexity  rapidly  increases  with  larger  J,  with 
frequency-selective  fading,  and  with  a  receiving  antenna  array  used  by  Bob  [5]. 


4.  Transmission  power  analysis 

Analysis  of  transmission  power  is  necessary,  not  only  for  enhancing  power  efficiency,  but 
also  for  guaranteeing  LPI.  Specifically,  we  need  both  to  reduce  the  total  transmission  power  and 
to  balance  the  power  among  the  transmitting  antennas.  We  will  show  in  this  section  that  this 
objective  can  be  conducted  by  choosing  properly  n  and  £ . 
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For  our  proposed  scheme,  conditioned  on  each  selected  channel  coefficient  ft, ,  the  total 
transmission  power  is 

tr{£[s(/?)sff  (/?)  |  h,P, (/?)]}  =  +  £}+|«,  |2  +f,H  +  £)f,  (1.9) 

whose  diagonal  entry  gives  the  transmission  power  of  each  antenna. 

Let  us  consider  specifically  the  case  that  fi  =  0  and  £  =  cr2IJ_1 .  The  total  transmission  power 
for  a  given  channel  realization  ft,  and  a  given  choice  of  ft,  becomes 

P,A  =  £'[si/(«)s(«)  |  h,P,(n)]  =  {J-\)cr2+  |  a,  |2  +||f,  ||2  a2  (1.10) 

Equation  (1.10)  shows  that  small  ft,  increases  the  total  transmission  power.  In  order  to  reduce 
transmission  power,  we  need  to  select  ft,  with  magnitude  larger  than  certain  threshold  a  ,  and  a 
should  be  carefully  selected.  Since  ft,  is  a  complex  Gaussian  random  variable  with  zero  mean 
and  unit  variance,  the  probability  for  the  selected  channel  coefficient  ft,  to  have  energy  greater 
than  a  is 

P[|  ft,  \2>  a]  =  \  e~'dt  =  e~a  (1.11) 

J  a 

In  other  words,  with  J  transmit  antennas,  the  average  number  of  selectable  coefficients 
is  Je~a  . 

Proposition  2.  With  Rayleigh  flat-fading  channels,  if  the  channel  coefficients  are  selected 
with  threshold  a  ,  then  the  expected  total  transmission  power  is 

Pt  =(/-l)o-2+l  +  (J-l)(l  +  o-2)r(0,«)  (1.12) 

From  (1.12),  we  can  see  that  the  expected  total  transmission  power  Pt  is  a  function  of  the 
number  of  transmitting  antennas  J ,  the  variance  <r2  of  the  random  weights,  and  the  threshold  a 
for  selecting  ft, .  Especially,  Pt  increases  when  a2  increases,  or  a  decreases,  or  J  increases. 
Figure  2  illustrates  their  relationships  under  /  =  4 . 


Figure  2.  Total  transmission  power  Pt  and  power  ratio  PuilPtJ  of  the  ith  transmit  antenna  to 
the  jtb  transmit  antenna  (j  *i)  when  ft,  is  selected.  Solid  lines:  total  power.  Dashed  lines: 

power  ratio. 
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If  the  channel  h  is  slowly  time-varying  or  even  constant  for  a  long  time,  we  need  to  avoid  the 
case  that  the  power  of  one  of  the  transmit  antennas  is  exceptionally  larger  than  the  others. 
Otherwise  the  array  behaves  as  a  single  antenna  and  the  security  can  be  compromised.  Therefore, 
we  have  to  constrain  the  ratio  of  the  transmission  power  of  the  ith  transmit  antenna 
Pt  i  =|  a,  |2  + 1|  f  ||2  a2  to  that  of  the  /"'  transmit  antenna  Ptj  =  a2 .  The  power  ratio  can  be  obtained  as 

Pu  l  +  (J-l)(l  +  cr2)r(0,a) 

P,.j 

Obviously,  it  is  usually  impossible  to  obtain  unit  ratio.  From  figure  2,  the  power  ratio  is  a 
decreasing  function  of  both  a2  and  a  . 

Larger  a2  increases  Pt  but  decreases  Pr,/PltJ  .  Since  both  Pt  and  PtJPtJ  should  be  small,  there 

is  a  trade-off  between  them  when  choosing  a2 .  In  the  simulations,  we  have  chosen  a2  =0.5.  On 
the  other  hand,  larger  a  reduces  both  Pt  and  Pu/PtJ  .  But  from  (1.11),  it  reduces  the  number  of 

selectable  ht  as  well  as  the  randomness  of  P,(«) .  Hence  there  is  also  a  trade-off  when  choosing  a  . 
We  have  used  a  =  0.5  in  simulations. 

5.  Simulations  and  experiments 

In  this  section,  we  use  three  simulation  experiments  to  study  the  effectiveness  of  the 
proposed  transmission  scheme  by  evaluating  the  BER  of  Bob  and  Eve.  Eve  is  assumed  to 
estimate  symbols  either  by  blind  equalization  (specifically,  via  the  CMA  algorithm  [11]),  or  by 
directly  using  Bob’s  method. 

In  the  first  simulation  experiment,  we  used  randomly  generated  channels.  For  comparison 
purpose,  we  evaluated  the  perfonnance  of  the  optimal  transmit  beamforming,  and  gave  the 
theoretical  BER  curve  of  the  Rayleigh  fading  channel  without  diversity.  The  simulation  results 
are  shown  in  figure  3.  Transmissions  with  the  proposed  Algorithm  1  have  similar  perfonnance  as 
the  optimal  transmit  beamforming.  Eve  can  not  intercept  symbols  using  blind  equalization.  On 
the  other  hand,  Algorithm  1  requires  both  larger  total  transmitting  power  and  larger  single¬ 
antenna  transmission  power  than  the  conventional  beamfonning.  In  addition,  the  standard 
deviation  of  transmission  power  becomes  large  as  well. 

Next,  considering  the  importance  of  verifying  the  extent  of  channel  similarity  between  Bob 
and  Eve,  in  the  second  simulation  experiment,  our  objective  is  to  show  how  confident  we  can 
claim  that  Bob  and  Eve's  channels  are  different.  We  considered  a  3x3x7  (height/wide/length,  in 
meters)  room  with  some  objects  (a  box  and  a  beam)  inside,  as  shown  in  figure  4.  We  placed  3 
transmitting  antennas  at  one  end,  and  523  receiving  antennas  at  the  other  end,  where  the 
receiving  antennas  were  put  on  a  grid  of  0.3  meter,  where  0.3  is  near  to  the  wavelength  of  1  GHz 
carrier.  Specifically,  there  were  15  planes,  each  had  35  receiving  antennas  (two  antennas  were 
missing  where  there  were  conflictions  with  the  objects).  We  let  the  transmit  antennas  to  transmit 
impulse  signals,  and  obtained  the  signals  received  by  each  of  the  receive  antennas.  This 
procedure  was  conducted  using  electromagnetic  (EM)  simulation  software  (based  on  FDTD). 
From  the  signals  we  then  estimated  all  effective  channels  on  a  0.3  meter  grid.  Details  of  the  EM 
simulation  and  source  data  can  be  obtained  at  [13]. 
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Figure  3.  Left:  Receiving  performance  comparison.  Right:  Total  transmission  power  and 
the  transmission  power  of  each  individual  antenna,  as  well  as  their  standard  deviations. 
Standard  deviation  is  shown  by  x  or  □  above  the  power  value. 


In  this  simulation,  we  obtained  altogether  523  array  channel  vectors  (J  =  3) .  Then  we  used 
each  of  them  as  Bob's  h  while  each  of  the  rest  as  Eve's  H(  to  examine  LPI.  Assuming  Eve  use 
Bob's  detection  method,  then  LPI  depends  on  the  difference  between  h  and  H .  Specifically,  the 
channel  difference  will  contribute  interference  to  Eve's  detection,  which  degrades  the  signal-to- 
interference  ratio  (SIR)  to  be  approximately 

SIR* - 7 - frf -  (1.14) 

|l  he  -h  ||2  yJT(J-l)cr2+ 1  a,  \2  + 1|  it  ||2  a2 
J  i= 1 

Therefore,  we  evaluated  the  cumulative  distribution  of  Eve's  SIR  (under  noiseless  assumption) 
for  523x522  possible  transmission/eavesdropping  cases.  The  results  are  shown  in  figure  4,  from 
which  we  clearly  see  that  in  almost  all  cases  (i.e.,  almost  100%),  Eve's  signals  suffer  a  very  high 
SIR  loss,  which  prevents  Eve  from  symbol  detection. 


Figure  4.  Left:  Settings  of  a  room  for  electromagnetic  wave  propagation  simulation.  Refer 
to  [13]  for  a  detailed  description.  Right:  Cumulative  distribution  function  of  the  SIR  of 
Eve's  signals  due  to  the  difference  between  Bob's  and  Eve's  channels.  Channels  are  derived 

by  EM  simulations. 
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As  the  third  experiment,  by  using  the  channels  obtained  from  the  EM  simulation,  we  have 
also  simulated  the  error  rates  of  Bob  and  Eve.  For  each  SNR  value,  Bob's  error  rate  was  the 
average  of  all  these  523  cases,  while  Eve's  error  rate  was  obtained  as  the  minimum  value  among 
all  523x522  cases  (100%)  or  the  majority  (99%)  of  523x522  cases.  The  results  are  shown  in 
figure  6  as  the  dashed  curves.  It  can  be  seen  that  for  almost  all  cases,  Eve's  error  rate  is  extremely 
large. 


Figure  5.  Experiment  setup  with  2  transmitting  antennas  (right)  and  2  receive  antennas 
(left).  Notice  the  short  distance  between  the  two  receive  antennas. 


O  Bob(EM) 
x  Eve(EM,100%) 
A  Eve(EM,99%) 
Bob(tcstbed) 
Eve(testbed) 


SNR  (dB) 


Figure  6.  BER  of  Bob  and  Eve.  Solid  lines:  using  channels  measured  from  testbed.  Dashed 
lines:  using  channels  obtained  from  EM  simulation. 


We  are  also  building  a  testbed  using  the  wireless  transmission  modules  of  ComBlock.com 
[13].  We  implemented  two  QPSK  transmitters  and  two  QPSK  receivers.  One  snap  shot  of  the 
experiment  is  shown  in  figure  5.  Four  channels  were  estimated  and  fed  into  the  program  of  the 
first  simulation  experiment  to  estimate  BER.  The  results  fit  well  with  those  obtained  by  purely 
simulations  (solid  lines  in  figure  6).  Note  that  the  two  receiving  antennas  (one  for  Bob,  one  for 
Eve)  were  purposely  placed  very  close  to  each  other. 
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Part  II:  Cooperative  Transmissions  with  Asynchronous  Transmitters 
1.  Introduction 

Recently  there  have  been  great  research  interests  in  cooperative  transmissions,  where  the 
increased  density  of  network  communication  devices  can  be  exploited  to  enhance  performance. 
Since  cooperative  transmissions  use  an  array  of  transmitters  similarly  as  physical  antenna  array 
transmissions,  the  rich  research  results  of  the  latter  can  naturally  be  extended  for  the  former. 
Especially,  space-time  block  codes  (STBC)  [14]  have  been  widely  investigated  for  cooperative 
transmissions  because  of  their  efficient  computation.  Nevertheless,  only  very  recently,  there  has 
been  work  toward  the  major  difference  between  cooperative  transmissions  and  conventional 
array  transmissions,  i.e.,  the  synchronization  of  cooperative  transmitters  [  1 5]-[  1 7] .  It  is  difficult, 
and  in  most  cases  impossible,  to  achieve  perfect  synchronization  among  distributed  transmitters. 

When  the  cooperative  transmitters  are  not  synchronized,  their  signals  arrive  at  the  receiver 
with  arbitrary  delays.  This  has  two  effects.  First,  any  intended  signal  structure,  such  as  the 
orthogonal  structure  of  STBC  encoding,  is  destroyed,  so  that  the  conventional  receivers  may  not 
work  anymore.  Next,  there  is  channel  dispersion  because  of  the  pulse  shaping  and  because 
universally  ideal  sampling  time  instants  for  all  transmitter's  signals  do  not  exist.  If  the 
environment  is  flat  fading,  as  expected  in  wireless  sensor  networks,  such  dispersion  makes  the 
channels  fundamentally  different.  In  contrast  to  flat  fading,  dealing  with  dispersive  channels  not 
only  makes  the  receiver  excessively  complex  and  energy  consuming,  but  also  greatly  reduces 
performance,  especially  if  the  optimal  maximum-likelihood  sequence  estimation  is  not 
affordable. 

Considering  that  it  is  still  not  very  clear  as  to  whether  cooperative  transmissions  really 
provide  positive  overall  gain  in  asynchronous  flat  fading  environment,  we  will  first  investigate 
the  performance  of  asynchronous  cooperative  transmissions  by  studying  the  ISI  due  to  channel 
dispersion,  and  by  comparing  its  signal-to-interference  ratio  (SIR)  against  the  original  flat  fading 
cases.  In  this  case,  we  mainly  consider  the  delay  asynchronism  and  ISI.  Next,  we  consider 
OFDM  transmissions  which  can  resolve  the  ISI  issues  efficiently,  but  suffer  from  the  carrier 
frequency  asynchronism. 


2.  Performance  of  cooperative  transmissions  with  delay  asynchronism  and  ISI 

We  consider  the  cooperative  transmission  among  nodes  1  to  J .  The  passband  signal  to  be 
transmitted  by  each  transmitter  /  has  a  general  form  Re| X, (w) PiM  ~ mT)ej2’lf  '] ,  where  si(m)  is 
the  complex  symbol  transmitted,  ph(t)  is  the  baseband  pulse  shaping  filter,  and  fc  is  the  carrier 
frequency.  After  delaying  with  5, ,  the  signal  received  by  the  receiver  is 


rp(t)  =  Re 


J  cc 

Tjai  X  sSm)Ph(t - mT - Si)e’1’lflAl-s,) 


i=\  m=- co 


+  v,(0 


(2.1) 


where  a,,  are  (complex)  fading  of  the  propagation.  Without  loss  of  generality,  we  can 
demodulate  (2.1)  to  obtain  the  continuous-time  complex  baseband  signal 
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(2.2) 


J  oo 

rb{t)  =  YjaieJ<>‘  Yj  s,(m)ph(t-mT -S^  +  v^t) 

1=1  m=-  oo 

where  the  phase  6i  =  -2k fc8i . 

Since  8,  may  be  different  for  different  transmitters,  without  loss  of  generality  we  perform 
baseband  sampling  at  time  instant nT  +  t  ,  which  gives  the  sample  x(n)  =  rh(nT  +  r) ,  where  0  <  r  <  T . 
The  sample  can  be  written  as 

./  co 

x{T\n)  =  ’^jaieie‘  ^  •s,.(w)t9i[(w-w)7’  +  r-£,]  +  v(r)(n)  (2.3) 

i= 1  /«-— co 

We  assume  that  the  pulse  shaping  fdter  have  support  [-KT,  KT] .  Obviously,  each  sample 
x(T\n)  depends  on  multiple  transmitted  symbols,  so  the  channel  should  be  dispersive. 

Define  channel  coefficients  (due  to  pulse  shaping  only,  not  the  propagation  flat  fading)  as 

hi{l)  =  Pb{lT  +  T-8l)  (2.4) 

Then  the  received  sample  can  be  written  as 

./  oo 

x(r,(/?)  =  ateJ0'  ^  hj(l)si(>i~l)  +  v(T)(n)  (2.5) 

1=1  l =-oo 

Define  the  lower  range  and  the  upper  range  of  /  as 


Then  the  channel  model  becomes 


(2.6) 


/  \s,(n-Lir) 

x(T\n)  =  YJaie]0'\Lhi{LLUr)  •••  A,(^r)]  :  +v(r)(«)  (2.7) 

Using  the  model  (2.7),  we  can  either  find  ways  to  mitigate  the  channel  dispersion,  and  hence 
ISI,  or  otherwise  suffer  from  ISI-induced  SIR  reduction.  To  evaluate  the  performance 
degradation  due  to  channel  dispersion,  we  consider  only  the  pre-processing  of  the  received 
samples  without  knowing  the  encoding  details.  We  will  try  to  restore  the  original  flat  fading 
channels  by  mitigating  dispersion  (though  we  still  have  delays  in  integer  symbol  intervals).  Then 
we  will  evaluate  the  interference  due  to  residue  dispersion,  and  the  reduction  of  SIR. 

As  a  basis  for  comparison,  we  let  the  receiver  just  use  the  flat  fading  channel  model  by 
assuming  all  other  channel  coefficients  in  (2.7)  are  interference.  We  define  the  SIR  as 


yJ  max,  |  /?.(/)  j2 

SIR  = - ,  ,  ^ -  (2.8) 

For  the  ideal  noisy  case  (without  asynchronism)  the  signal-to-noise  ratio  (SNR)  can  be 
defined  as  SNR  =  Jo]  jo] ,  which  is  also  the  average  over  the  fading  coefficients.  Then,  in  case  of 


asynchronism,  the  extra  interference  in  (2.7)  makes  the  receiver's  signal-to-interference-and- 
noise  ratio  (SINR)  to  be 
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SINR  = 


1 


1 


1 


(2.9) 


./SIR 


SNR 


Therefore,  both  SIR  and  the  reduced  value  of  channel  coefficients  max,  |  ht(l)  |2  make  the 
SINR  lower  than  the  ideal  case  SNR.  In  particular,  if  SIR«SNR,  then  SIR  dominates  and 
becomes  the  SINR  floor,  which  means  in  high  SNR  environment,  the  SINR  is  limited  to  be  SIR. 
On  the  other  hand,  if  SIR»SNR,  then  the  interference  does  not  degrade  system  performance  too 
much. 


By  examining  SIR  instead  of  SINR,  we  can  reduce  complexity  while  are  still  able  to  show 
the  impact  of  the  asynchronism-induced  interference.  Therefore,  to  simplify  the  problem,  we 
consider  noiseless  transmission  so  as  to  focus  on  the  interference. 


2.1  Use  fixed  sampling  timing 

We  refer  (2.7)  as  the  “fixed”  sampling  scheme,  and  use  it  as  a  basis  for  comparison.  With  the 
fixed  scheme,  without  loss  of  generality,  we  can  use  r  =  0 .  We  can  decompose  the  delay  8,  into 
an  integer  value  Sf  and  a  fractional  value  S'f .  In  addition,  we  assume  that  the  delays  of  the 
transmitters  are  independent  and  uniformly  distributed,  which  means  S'f  is  a  uniform  random 
variable.  Then  SIR  (2.8)  becomes 


SIR  = 


yIwlM(*r+3(2)n  I2 

where  the  range  of  the  integer  k  in  the  denominator  of  (2.10)  is 

\-K-8f~\<k<[_K-5f\,  k*  0 


(2.10) 


(2.11) 


2.2  Mitigate  interference  by  optimizing  sampling  timing 

A  simple  enhancement  scheme  is  to  look  for  the  optimal  time  instant  r„  such  that 
r0  =  arg  max  SIR.  With  the  optimized  ro  ,  we  have  channel  coefficients  h,.(l)  =  pb{lT  +  z0  -8t)  . 

re[0,r) 

Similarly,  we  can  define  r0  -8j  =  dT  +  zT  ,  where  d,  is  an  integer  andz,  e  (-1/2, 1/2] .  Then  the  SIR 
expression  (2.8)  becomes 


SIR  = 


ZiljMAOl 


yZwZM(*T  +  z/T)  I2 

where  the  range  of  the  integer  k  in  the  denominator  of  (2.12)  is 


(2.12) 


\-K-zf<k<lK-Zi],  k*  0  (2.13) 

The  major  difference  between  (2.10)  and  (2.12)  is  that  8'f'  has  uniform  distribution,  but  z, 
has  a  distribution  more  concentrated  around  0,  so  that  (2.12)  is  more  likely  larger  than  (2.13). 


2.3  Mitigate  interference  by  over-sampling  and  combining 
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The  basic  idea  is  that  the  received  signals  have  a  rich  structure  because  of  different  delays,  so 
the  over-sampled  signals  contain  new  infonnation.  Consider  over-sampling  by  a  factors,  which 
means  we  have  iV  +  1  samples  in  each  symbol  interval  at  time  instants  r(: ,  k  =  0,1, •••,7V ,  where 
0  <  rk  <  T .  Then  for  each  rk ,  we  have  a  sample  model  similar  to  (2.7).  We  can  stack  JV  +  1  samples 
of  each  symbol  interval  for 


x(«)  =  a,eJV<  H,s,  (/?)  +  v(«) 

1=1 

where  the  (N  + 1)  x  (If  -L'  +1)  dimensional  channel  matrices 


H  = 


,  WL) 


0 


h,(LlK)  -  A,(i£„)  0if_, 


Note  that  0,  denotes  a  lx£  dimensional  zero  vector. 


(2.14) 


(2.15) 


Then  we  can  combine  the  samples  by  an  N  + 1  dimensional  vector  f  to  obtain  y(n)  =  f"x(») . 
Ideally  (with  zero-forcing  criterion)  we  expect 

y(n)«Y,aiej0‘gisi(n-di)  +  win)  (2.16) 

i=l 

which  becomes  a  flat  fading  channel  model,  although  there  may  have  integer  delays.  The  key 
point  is  that  the  problems  with  respect  to  dispersive  channels  are  gone. 

Using  zero-forcing  criterion  means  thatf"H  /  =  1, ■■•../ ,  where  erf  is  a  unit  vector  with 

value  1  in  the  d''1  entry  and  zeros  elsewhere.  One  way  to  calculate  the  vector  f  is 


f"=[<  -  <][H,  ...  H,]*  (2.17) 

The  SIR  of  y(n)  can  be  calculated  similarly  as  (2.8),  but  with  the  composite  vectors  f"H,  in 
place  of  h,.(l) . 


2.4  Simulations  and  performance  comparison 

In  this  section,  we  use  simulations  to  compare  the  SIR  of  the  three  schemes:  fixed  sampling 
time  (Section  2.1),  selective  optimized  sampling  time  (Section  2.2),  and  linear  combining 
(Section  2.3).  Figure  7  shows  the  cumulative  distribution  of  SIR.  Specifically,  we  see  that  for 
fixed  method,  the  SIR  has  80%  probability  of  being  lower  than  15  dB.  But  for  the  selective 
method,  the  SIR  has  80%  probability  to  be  higher  than  15  dB.  The  combining  method  has  99% 
probability  higher  than  30  dB. 
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Figure  7.  Cumulative  distribution  of  SIR. 


3.  Using  OFDM  transmissions  to  combat  channel  dispersion 

Orthogonal  frequency  division  multiplexing  (OFDM)  transmission  technique  is  desirable  for 
combating  the  loss  of  timing-phase  synchronization,  since  any  limited  propagation  delay  (or 
timing-phase)  difference  among  the  signals  of  cooperative  transmitters  can  be  tolerated  by 
simply  increasing  the  length  of  cyclic  prefix  (CP)  [18]  [19].  Because  of  this,  it  may  find  wide 
applications  in  cooperative  transmissions.  Nevertheless,  OFDM  suffers  critically  from  the  loss  of 
carrier  frequency  synchronization  where  the  carrier-frequency  offset  (CFO)  incurs  inter-carrier 
interference  (ICI).  This  CFO  problem  becomes  much  worse  in  multi-transmitter  OFDM  systems 
because  of  the  increase  in  inter-transmitter  interference,  not  only  ICI  [19].  Unfortunately,  there 
have  been  no  method  that  claims  complete  CFO  cancellation  (to  our  knowledge),  and  most 
existing  methods  primarily  deal  with  small  CFO  only. 

Considering  that  the  CFO,  or  the  overhead  of  achieving  perfect  synchronization  otherwise, 
may  severely  degrade  the  gain  of  cooperation,  it  is  more  desirable  to  develop  methods  for  the 
receiver  to  cancel  CFO  completely  in  cooperative  OFDM  transmissions.  We  here  present  a  novel 
approach  to  realize  such  an  objective.  Our  basic  idea  is  to  utilize  the  redundancy  of  the  long  CP 
for  CFO  mitigation  or  even  cancellation.  A  unique  feature  of  our  approach  is  that  it  is 
implemented  purely  as  a  “pre-processing”  procedure,  independent  from  cooperative 
encoding/decoding  details.  In  other  words,  it  simply  makes  the  CFO  problem  transparent  to  the 
cooperative  OFDM  transmission  designs. 

Consider  a  cooperative  transmission  system  with  I  cooperative  transmitters  and  one  receiver. 
All  the  I  cooperative  transmitters  are  assumed  to  have  the  same  data  packet  that  is  to  be  encoded 
and  transmitted,  using  some  predefined  cooperative  encoding  schemes  such  as  cooperative 
STBC.  The  encoder  outputs  Z>.(«)  are  then  OFDM  modulated,  which  gives  the  OFDM  signal  st(n) . 
The  discrete  baseband  channel  from  the  /th  transmitter  to  the  receiver  is  assumed  frequency 
selective  fading  with  coefficients  h,(l) . 
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Figure  8.  Multi-transmitter  cooperative  OFDM  transmission  and  receiving  block  diagram. 

From  the  received  signal  r(n) ,  the  receiver  mitigates  the  asynchronism  in  carrier  frequency 
and  timing  using  our  proposed  method,  after  which  conventional  OFDM  demodulation  and 
cooperative  decoding  techniques  are  applied. 

With  the  consideration  of  asynchronous  transmitters,  the  signal  of  each  transmitter  i  may 
have  a  propagation  delay  ds  and  a  CFO  si  (relative  to  a  reference  timing  and  a  reference  local 
carrier)  when  received  at  the  receiver.  We  assume  di  to  be  integer  (with  symbol  interval  as  unit) 
since  the  fractional  portion  of  the  delay  contributes  nothing  but  some  extra  channel  dispersion 
which  can  be  assimilated  into  the  dispersive  channel  model.  The  CFO  ei  is  derived  as  the 
residual  carrier  frequency  normalized  by  the  OFDM  sub-carrier  frequency  separation. 

The  transmitted  signal  s,(n)  is  derived  from  the  Inverse  Fourier  Transform  (IFT)  of  the 
encoded  symbol  /?,.(«)  •  Since  there  is  no  inter-block  interference  (IBI)  thanks  to  cyclic  prefix,  we 
consider  one  OFDM  block  for  notational  simplicity.  Then  the  /th  transmitter's  signal  st(n)  can  be 
written  as 


si(n)  =  -}=Yjbi(k)ei2*nk/N,  -Ng<n<N- 1  (3.1) 

k= 0 

where  Ng  is  the  length  of  the  CP  and  N  is  the  IFT  block  length  (we  also  define  it  as  OFDM 
block  length).  The  composite  signal  received  by  the  receiver,  with  delay  d,  and  CFO  ei 
considered,  is 

r(0)  =  ]TeME,.(0)H|.sl.  +  v(0)  (3.2) 

1=0 

where  $  is  the  initial  phase,  the  sample  vector  r(0)  =  [r(0)  •••  r(7V  —  l)]r ,  the  NxN  diagonal 
matrix 

E;. (0)  =  diag  jl, ejs‘ eJe‘ (AM) j  is  the  CFO  matrix,  the  symbol  vector  s,.  =  [si(0),---,si(N-l)f ,  and  the 
channel  matrix  is  NxN  circulant 


~  0W_W, 

ML) 

■■  MO) 

°4 

«,= 

"N-L-d,+ 1 

ML)  - 

■■  MO) 

0, 

^N-L-d,- 1 

ML)  ■ 

■■  MO) 

0, 

Our  basic  idea  is  to  exploit  the  redundancy  of  the  CP  based  on  the  structure  of  the  signal 
model  (3.2).  If  the  CP  length  N  is  longer  than  L  +  max  d  ,  then  in  addition  to  those  inr(0),  we 

g  0<i</-l 

have  more  IBI-free  samples  with  which  we  can  construct  new  sample 
vectors  r (m)  =  [r(-m), •  •  ■ ,  r(N  - 1  -  m)f .  Similarly  to  (3 .2),  we  have 
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(3.4) 


where 


/-I 

r  (m)  =  ^V'^E^mjFTs,  +  v(/w) 

1=0 


E  ,('«)  = 


0 


(w|Af)x(Af-m|iV) 


diag{e 


JS,  (-'») 


•,e 


diag{e 


j£j  ( -m+m\N ) 


•,e 


0 


(iV-m|Af)x(m|Af) 


Note  that  (x  |  At)  denotes  x  mod  N . 


(3.5) 


Noticing  that  (3.2)  and  (3.4)  contain  the  same  H  and  s,  but  have  different  CFO  matrices,  we 
can  stacking  together  all  available  vectors  r(m) , 


r(0) 

7-1 

~E,(0)~ 

H/S ,  + 

v(0) 

r(M-l) 

i=0 

E,(0)_ 

v(M-l) 

which  for  notational  simplicity  can  be  denoted  as 

/-I 

y  =  ^e^'A.H.s,.  +u 

i=0 

The  dimensions  of  y  and  A,  are  MNx l  and MN x  ,v ,  respectively. 


(3.6) 


(3.7) 


Our  basic  idea  is  thus  to  design  an  NxMN  CFO  mitigation  matrix  X  such  that 

xa,=iv  (3.8) 

for  all  i  =  0, •■■,1-1.  If  X  is  available  for  (3.8),  then  CFO  can  be  mitigated  via 


z  =  Xy 

Note  that  a  straightforward  solution  for  X  is 


(3.9) 


x  =  [iff 


l»] 


E„(0) 

E0(M-1) 


E,_,  (0) 


E;_j  (M-l) 


(3.10) 


If  (3.8)  can  be  satisfied  perfectly,  then  we  havez  =  =1e^iHjsi  +Xu ,  which  is  a  conventional 

CFO-free  OFDM  sample  vector  after  removing  the  CP.  With  the  vector  z ,  conventional  OFDM 
demodulation  can  be  applied  to  detect  symbols  b,(k) . 


We  have  found  in  [P2]  that  CFO  can  be  completely  removed  only  if  CP  is  long  enough  and 
appropriate  sample  vectors  r(m)  are  used.  Not  all  available  r (»;)  need  to  be  used  though  and  in 
fact,  using  less  r(m )  leads  to  reduced  complexity.  The  inverse  of  the  big  matrix  X  can  be 
avoided  by  exploiting  the  special  structure  ofE,(m) . 

Specifically,  the  computational  complexity  consists  of  two  parts.  The  first  part  is  the 
calculation  of  X ,  where  the  good  news  is  that  X  needs  to  be  calculated  only  once  (for  all  OFDM 
blocks)  if  si  is  not  time-varying.  In  this  case,  the  complexity  is  in  the  order  of  0(NI3) .  The  second 
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part  is  the  calculation  of  Xy  for  each  OFDM  block,  where  the  complexity  is  0(  NI)  since  the 
majority  of  X  entries  are  zeros.  Considering  that  the  first  part  happens  only  once  and  /  (the 
number  of  cooperative  transmitters)  is  usually  much  smaller  than  .V ,  the  proposed  algorithm  has 
complexity  almost  linear  in  N ,  which  is  very  efficient. 

We  have  also  derived  the  following  proposition  regarding  complete  CFO  cancellation. 

Proposition  2:  Let  the  CP  length  beN  >  ( /  —  1 ) N  +  L  +  max  di .  With  appropriate  parameters  m 

8  0<i<I 

forr(m) ,  CFO  can  be  completely  canceled  if  ei  -ej  *  InljN ,  for  any  i  *  j  and  integer  / . 

The  proposed  method  is  desirable  for  cooperative  transmissions  in  ad  hoc  wireless  networks, 
where  the  long  CP  (repeated  transmissions  like  spectrum-spreading)  is  used  for  CFO  cancellation, 
for  high  transmission  power  efficiency  as  well  as  for  better  noise/interference  suppression.  Some 
other  benefits  of  the  implementation  specified  by  the  Proposition  2  come  from  the  Vandermonde 
matrix.  Vandennonde  equation  systems  have  very  efficient  algorithms  to  solve,  with  complexity 
0(I2)  instead  ofO(/3).  As  a  result,  the  complexity  of  calculating  X  becomes  0{NI2) ,  instead 
of  0(NI3)  .  Furthermore,  Vandermonde  system  solver  can  usually  give  surprisingly  accurate 
solutions,  even  for  ill-conditioned  matrix.  This  property  is  especially  helpful  in  the  case  where 
some  CFOs  si  are  close  to  each  other. 

In  order  to  evaluate  the  performance  of  our  algorithm,  we  have  simulated  a  system  with  two 
cooperative  transmitters  and  one  receiver,  using  Alamouti  STBC.  We  compared  the  perfonnance 
of  our  algorithm  against  the  ideal  cooperative  transmissions  with  perfect  synchronization 
(“Perfect”),  as  well  as  two  other  OFDMA  CFO  mitigation  methods:  [20]  (“CLJL”)  and  [21] 
(“HL”).  We  used  N  =  32 ,  QPSK.  The  integer  delays  di ,  the  CFOs  st ,  and  the  channels  (with 
order  L  =  3)  were  all  randomly  generated  for  each  transmitter  during  each  run  of  the  simulation. 
We  used  10,000  runs  of  the  simulations  to  derive  the  average  symbol  error  rate  (SER)  under 
various  signal-to-noise  ratios  (SNR)  or  various  CFO. 

From  the  simulations  results  shown  in  figure  9,  we  can  see  that  our  algorithm  has  much 
better  performance  when  SNR  is  not  extremely  low,  and  the  advantage  is  even  more  significant 
for  large  relative  CFO.  In  particular,  when  relative  CFO  is  0.5,  “HL”  and  “CLJL”  failed,  but  our 
algorithm  had  a  performance  almost  independent  of  the  size  of  relative  CFO.  In  addition,  we 
have  verified  the  “independence”  of  our  algorithm  to  various  CFOs.  The  conventional  method 
did  not  resolve  the  CFO  problem,  and  neither  did  the  “HL”  scheme  when  the  rCFO  was  not  very 
small.  In  contrast,  our  new  method  showed  almost  constant  perfonnance  under  various  rCFO. 
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Figure  9.  Left:  Performance  comparison  of  our  algorithm  with  HL  [21]  and  CLJL  [20]  for 
CFO  mitigation  under  rCFO  0.1  and  0.5.  Right:  Performance  comparison  of  our  “New” 
CFO  mitigation  algorithm  with  the  conventional  OFDM  receiver  and  the  CFO  mitigation 

algorithm  “HL”  [21].  SNR  20  dB. 

4.  Extensions  to  MC-DS-CDMA 

In  [P3],  we  have  extended  the  CFO  mitigation  method  in  Section  3  into  MC-DS-CDMA 
systems,  and  have  proposed  a  new  receiving  algorithm  for  MC-DS-CDMA  systems  when  carrier 
frequency  offset  (CFO)  is  significant.  By  exploiting  the  special  structure  of  the  CFO 
contaminated  signals,  the  new  algorithm  cancels  CFO  completely  during  the  dispreading 
procedure,  after  which  the  despread  CFO-free  signal  is  demodulated  via  normal  FFT-based 
OFDM  demodulator.  While  guaranteeing  complete  CFO  cancellation,  this  method  is 
advantageous  over  the  majority  existing  CFO-mitigation  techniques  that  can  only  mitigate  but 
not  completely  remove  CFO.  An  efficient  algorithm  is  developed,  and  simulations  are  conducted 
to  demonstrate  the  performance. 

Cooperative  communication  is  an  important  concept  in  wireless  networks.  We  will  continue 
study  its  implementations  and  applications,  not  only  in  wireless  security,  but  also  in 
communication  efficiency.  On  the  other  hand,  we  have  also  noticed  its  interesting  relation  to 
another  important  concept:  cross-layer  design.  Work  is  being  conducted  by  us  along  this  line. 
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Part  III:  Cooperative  Transmissions  Testbed  Development 


1.  Introduction 

The  objective  for  this  testbed  development  is  to  demonstrate  the  concepts  of  cooperative 
transmissions  and  secure  array  transmissions  that  we  have  developed.  By  using  wireless 
communication  modules  purchased  from  comblock.com,  we  are  building  a  wireless  transmission 
testbed  for  demonstrating  our  results  in  both  secure  communications  and  cooperative 
communications.  As  specified  in  the  research  topic  4.2  of  the  project  proposal,  the  testbed  is 
needed  to  show  that  the  proposed  security  schemes  work  in  practice  and  also  to  show  how 
difficult  it  is  to  implement  it.  The  testbed  is  also  important  to  demonstrate  the  cooperative 
transmission  schemes,  since  we  have  considered  practical  problems  such  as  the  lack  of 
synchronization. 

The  reason  that  we  use  comblock.com  modules  is  for  low-cost  implementation.  However, 
those  devices  are  for  single  antenna  transmission  and  receiving  only,  not  for  antenna  array 
transmissions.  Therefore,  we  need  overcome  this  major  hurdle,  i.e.,  find  ways  to  use  such 
ComBlock  modules  to  emulate  antenna  array  transmissions  and  receiving.  We  have  been  trying 
two  different  approaches  for  this  objective. 

The  first  approach  is  to  re- wire  hardware.  After  some  studying  and  trying,  we  find  a  specific 
wire  that  initializes  the  working  of  the  module.  Therefore,  we  replace  this  wire  by  an  external 
trigger  signal  that  is  under  our  control,  by  which  we  are  able  to  initiate  the  transmission  of  the 
two  transmitters  at  the  same  time,  so  as  to  have  them  to  emulate  antenna  array  transmissions. 

The  second  approach  is  to  develop  new  software  to  replace  the  existing  “ComBlock  Control 
Center”  software  provided  by  comblock.com.  The  existing  software  is  designed  to  control  one 
ComBlock  module  only,  so  for  two  transmitters,  we  have  to  use  either  two  computers  or  to  open 
two  “ComBlock  control  center”  windows  in  the  same  computer.  This  means  it  is  impossible  for 
us  to  make  the  two  transmitters  to  begin  transmit  at  the  same  time.  Our  objective  is  to  design 
new  software  so  that  we  can  initiate  the  two  transmitters  at  the  same  time.  Much  effort  has  been 
spent  from  the  summer  of  2006,  and  some  initial  progress  has  been  achieved.  For  example,  we 
have  setup  the  new  control  window,  and  can  “ping”  the  modules  via  TCP/IP  communications. 
We  have  also  demonstrated  that  our  software  can  initiate  the  transmission  of  two  transmitters, 
which  is  achieved  by  writing  to  some  special  registers.  Nevertheless,  this  software  approach 
needs  more  time  to  finish.  So  far,  most  of  our  efforts  are  in  the  hardware  approach. 

We  have  adopted  a  hardware  approach  to  wire  the  two  transmitters  together.  In  order  for  the 
two  transmitters  to  transmit  at  approximately  the  same  time,  we  use  the  optional  external  trigger 
signal  of  the  ComBlock  devices,  and  use  an  external  switch  to  initiate  the  transmission  of  the  two 
devices  at  the  same  time.  In  addition,  the  two  transmitters  need  to  guarantee  a  similar  carrier 
frequency,  which  is  achieved  by  using  a  common  external  carrier  frequency  reference  signal. 
Another  issue  is  the  sampling  clock,  which  fortunately  does  not  show  any  extra  requirement  for 
further  synchronization.  The  sampling  clocks  of  the  two  transmitter  devices  are  already  close 
enough.  But  if  otherwise,  we  can  also  use  an  external  sampling  clock  reference  signal.  So  far  this 
does  not  seem  necessary. 
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2.  Testbed  development  description 

First,  we  have  screwed  the  devices  onto  hard  plastic  boards.  Each  whole  transmitter  or 
receiver  is  screwed  on  one  board,  as  shown  in  figure  10. 


Figure  10.  Left:  two  transmitters.  Right:  two  receivers. 

Next,  in  order  to  explain  the  working  of  the  transmitters,  let  us  first  have  a  look  at  the  whole 
transmitter  structure  shown  in  figure  1 1 . 


Using  the  ComBlock  Control  Center  GUI  software  provided  by  comblock.com,  the 
initialization  of  the  transmission  is  as  follows:  1)  upload  the  transmission  data  (generated  from 
MATLAB  modulator)  into  Com8001;  2)  in  ComBlock  Control  Center,  select  “continuous 
download”  to  begin  transmission. 

We  exploit  the  fact  that  ComBlock  devices  provide  an 
optional  initiation  approach,  i.e.,  using  an  external  trigger  signal. 

In  Com8001,  there  is  an  “EXTTRIGGERIN”  input  signal 
(connected  to  Com5001  EXT  TRIGGER  OUT).  This  trigger 
signal  is  optional  in  default.  To  use  this  external  trigger  feature, 
we  can  apply  either  a  software  approach  or  a  hardware  approach. 

We  have  thus  far  implemented  the  hardware  approach.  As 
shown  in  the  right  figure,  we  connect  the  EXT  TRIGGER  IN 
to  a  3.3V  output  signal  provided  by  the  Com8001,  through  a 
switch.  So,  in  order  to  initiate  transmission,  we  do  the  following: 

1)  Upload  the  transmission  data  (generated  from  MATLAB 
modulator)  into  Com8001; 
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2)  Use  ComBlock  Control  Center  to  set  Com8001  in  “continuous  download”  and  “external 
trigger  enable”  modes,  i.e.,  set  the  Register  19  to  “09”  for  the  former  and  the  Register  21 
to  “01”  for  the  latter; 

3)  Press  the  switch  once,  then  the  transmitter  begins  transmission  repeatedly.  Note  that  in 
order  to  stop  the  transmission,  we  just  need  to  use  ComBlock  Control  Center  to  set  the 
Register  19  of  Com8001  into  “OB”. 

The  same  trigger  signal  can  be  sent  to  anther  transmitter,  so  that  the  two  transmitters  can  be 
set  to  start  transmission  at  approximately  the  same  time.  This  is  the  hardware  approach  for  us  to 
wire  the  two  independent  transmitters  for  emulating  a  transmission  array. 

In  order  to  make  sure  that  the  two  transmitters  have  approximately  the  same  carrier 
frequency,  we  can  use  the  same  external  10  MHz  reference  signal.  The  Com4001  can  use  either 
the  internal  reference  signal  or  the  external  10  MHz  reference  signal  for  the  oscillator  to  generate 
the  carrier  signal.  As  shown  in  figure  1 1,  we  provide  the  external  10  MHz  reference  signal  (from 
a  signal  generator)  to  the  Com4001  of  both  transmitters.  Then  we  use  ComBlock  Control  Center 
to  set  the  Com4001  Register  6,  bit  1,  to  be  “1”  to  use  the  external  reference  signal.  Note  that  this 
is  for  emulating  antenna  array  transmissions.  For  emulating  cooperative  transmissions,  we  may 
not  need  such  a  common  reference  signal,  because  carrier  frequency  difference  is  actually 
something  we  have  to  deal  with,  as  shown  by  our  research  results  described  in  the  last  report. 

With  the  external  reference  signal,  the  carrier  frequencies  of  the  two  transmitters  become 
very  similar.  In  our  preliminary  modulation/demodulation  experiments,  the  two  transmitters  have 
carrier  frequencies  as  close  as  to  100  Hz.  Furthermore,  this  small  discrepancy  may  be  due  to  our 
current  demodulation  algorithm  that  does  not  estimate  carrier  frequency  accurately  enough. 

Using  the  above  wiring  schemes,  we  have  also  tried  array  transmissions.  The  demodulation 
results  are  very  good,  almost  without  errors.  The  error  rate  is  only  slightly  worse  than  a  single 
transmitter  case.  Note  that  this  preliminary  result  only  shows  that  the  array  transmission  is  in 
effective.  In  order  to  estimate  channels  from  two  transmitters  by  only  a  single  receiver,  we 
simulated  direct-sequence  spread  spectrum  (DSSS)  transmissions,  and  compiled  receiving 
algorithms.  We  have  successfully  run  the  experiments,  and  have  estimated  the  channels.  Detailed 
of  experiment,  MATLAB  codes,  and  source  data  can  be  found  in  [13]. 

3.  Some  observations  and  future  work 

The  testbed  is  very  good  to  demonstrate  the  cooperative  transmissions  with  asynchronous 
transmitters. 

We  have  clearly  seen  the  delay  difference  and  carrier  frequency  difference  between  the  two 
transmitters.  We  have  two  journal  papers  in  preparation,  one  for  the  synchronization-induced  ISI, 
and  the  other  for  CFO  mitigation  in  MC-DS-CDMA.  We  have  the  plan  of  including  the  testbed 
experiments  as  demonstrations  in  these  two  papers. 

In  addition,  we  have  found  that  the  residue  phase  of  the  carrier  has  a  significant  impact  on  the 
channels.  We  have  run  multiple  experiments,  and  found  that  the  channels  usually  have  the  same 
power,  but  different  phase  from  each  running.  This  indicates  that  propagation  channels  are  quite 
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constant,  but  they  have  a  phase  changing  due  to  the  phase  synchronization  error.  This  effect  has 
both  advantages  and  disadvantages  for  our  proposed  secure  transmissions. 

The  advantage  is  that  the  phase  difference  can  be  exploited  by  the  transmitters  and  the 
receivers  to  make  the  channels  even  more  independent.  For  example,  some  intentional  phase 
jitter  may  be  introduced  by  the  PLL  circuits  to  make  the  channels  different. 

On  the  other  hand,  the  disadvantage  is  that  a  time -varying  channel  may  change  the  channel 
reciprocity,  or  may  prevent  the  application  of  channel-feedback  for  the  transmitters  to  get  and 
utilize  the  channel  state  infonnation. 

We  have  realized  that  this  special  phase  synchronization  problem  is  specific  to  our  testbed. 
Because  the  receiver  has  no  built-in  PLL  for  carrier  tracking,  but  just  using  its  own  local  carrier 
to  conduct  demodulation,  it  relies  on  the  receiving  program  to  compensate  for  the  carrier  offset. 
We  have  programmed  successfully  to  track  perfectly  the  carrier  frequency,  but  the  phase  tracking 
is  still  a  problem. 

In  the  future  work,  we  will  continue  to  work  on  the  phase  tracking  problem,  and  would  like 
to  see  how  accurate  we  can  track  the  phase.  Although  we  have  found  that  our  program  can 
successfully  track  phase  in  clear  (good)  artificial  signals,  but  it  still  has  some  problems  when 
applied  to  real  received  signals.  Furthermore,  we  would  like  to  see  how  actual  analog/digital  PLL 
circuits  such  as  the  Costas’  Loop  can  track  the  phase. 

After  we  have  the  new  phase  tracking  programs  that  are  able  to  process  real  signals,  we  will 
estimate  the  channel  properties  of  the  array,  and  test  our  secure  transmission  schemes.  In 
particular,  we  need  to  verify  the  channel  reciprocity  or  channel  feedback  property.  Various 
cooperative  transmission  schemes  can  also  be  verified. 

A  software-approach  can  be  helpful  for  more  flexible  experiments,  not  only  for  our 
secure/cooperative  transmissions,  but  also  for  extending  our  testbed  to  a  more  general  software- 
defined  radio  or  cognitive  radio.  For  this  purpose  we  need  to  program  our  own  control  software 
instead  of  using  the  ComBlock  Control  Center.  So  far,  we  have  showed  the  critical  TCP/IP 
communications  between  our  program  and  the  ComBlock  devices.  As  a  result,  the  subsequent 
work  might  just  be  some  painstaking  programming. 
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List  of  Acronyms 


AFRL 

Air  Force  Research  Laboratory 

AGWN 

Additive  Gaussian  White  Noise 

BER 

Bit  Error  Rate 

BPSK 

Binary  Phase-Shift  Keying 

CDMA 

Code  Division  Multiple  Access 

CFO 

Carrier  Frequency  Offset 

CP 

Cyclic  Prefix 

DSSS 

Direct  Sequence  Spread  Spectrum 

EM 

Electromagnetic 

FDTD 

Finite  Difference  Time  Domain 

FFT 

Fast  Fourier  Transform 

GUI 

Graphical  User  Interface 

IBI 

Inter-Block  Interference 

ICI 

Inter-Carrier  Interference 

IFT 

Inverse  Fourier  Transform 

LPI 

Low  Probability  of  Intercept 

MILCOM 

Military  Communications  Conference 

OFDM 

Orthogonal  Frequency  Division  Multiplexing 

OFDMA 

Orthogonal  Frequency-Division  Multiple  Access 

OPSK 

Quadrature  Phase-Shift  Keying 

PI 

Principal  Investigator 

PLL 

Phase  Locked  Loop 

SER 

Symbol  Error  Rate 

SINR 

Signal-to-Interference-and-Noise  Ratio 

SIR 

Signal-to-Interference  Ratio 

SNR 

Signal-to-Noise  Ratio 

STBC 

Space-Time  Block  Codes 

TCP/IP 

Transmission  Control  Protocol/Internet  Protocol 
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